Spin-Depairing Transition of Attractive Fermi Gases on a Ring 
Driven by Synthetic Gauge Fields 
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Motivated by the recent experimental realization of synthetic gauge fields in ultracold atoms, we 
investigate one-dimensional attractive Fermi gases with a time-dependent gauge flux on the spin 
sector. By combining the methods of the Bethe ansatz with complex twists and Landau-Dykhne, it 
is shown that a spin-depairing transition occurs, which may represent a nonequilibrium transition 
from fermionic superfluids to normal states with spin currents caused by a many-body quantum 
tunneling. For the case of the Hubbard ring at half filling, our finding forms a dual concept with 
the dielectric breakdown of the Mott insulator discussed in Phys. Rev. B 81, 033103 (2010). We 
analyze cases of arbitrary filling and continuum model, and show how filling affects the transition 
probability. 
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I. INTRODUCTION 

Owing to high controllability and flexibility of exper- 
iments, ultracold atoms have offered a testing ground 
to simulate a strongly correlated system in one dimen- 
sion both under equilibrium and nonequilibrium condi- 
tions Jj, ^2,] . Recent experimental efforts make it possible 
to study the system on a ring, which has been realized 
in continuum space [1, [3] and in optical lattice @ . As 
far as two-component Fermi gases on a ring in equilib- 
rium are concerned, one may obtain exact solutions by 
the Bethe ansatz method regardless of continuum Q or 
lattice 0, and many of the properties of the systems 
can be examined in combination with the bosonization 
approach. However, as for what happens in the Fermi 
gases in nonequilibrium, a unified approach has yet to 
be established and our current knowledge is far from a 
comprehensive level. 

Meanwhile, in recent years studies of synthetic gauge 
fields have attracted attention in cold-atom community 
P and the NISTgroup has succeeded in producing syn- 
thetic magnetic Q and electric fields 10], and spin-orbit 
coupling 'll'] on a continuum space. In particular, the 
electric fields discussed in Ref. [Ifl] have been realized 
with the creation of time-dependent gauge fields. In addi- 
tion, very recently, Aidelsburger et al. have succeeded in 
creating synthetic gauge fields on an optical lattice [l^ . 
which may pave the way to areas of study in cold atomic 
gases. In this paper, based on the situations described 
above, we consider a two-component attractive Fermi gas 
on a ring with a time-dependent gauge fiux on the spin 
sector at absolute zero. In the absence or presence of a 
time-independent gauge flux, it is well known that the 
ground state of the two-component attractive Fermi gas 
is filled with bound states of up-spin and down-spin parti- 
cles and the spin excitation has a gap, which is attributed 
to the appearance of fermionic superfluidity [is!]. Here, 
we demonstrate that in the presence of a time-dependent 
gauge flux <I>(t) = Ft, a spin-depairing transition induced 
by a many-body quantum tunneling occurs, which may 



represent a breaking phenomenon of fermionic superfluid- 
ity in nonequilibrium. Our prediction relies on a recently 
proposed method [ij, [l^ combining the Bethe ansatz 
method with complex twists [l^ [iJl and the Landau- 
Dykhne method |18l420l [ with which a quantum tunneling 
probability between two states can be evaluated. In fact, 
for the Hubbard model, our flnding builds a dual concept 
with the nonequilibrium transition from Mott insulators 
to metals in the repulsive Hubbard ring at half flUing 
with a time-dependent charge flux discussed in the con- 
text of electronic systems IJ] . In addition, we show that 
by using the formalism with dressed energies [l^ [2l| . 
the spin-depairing transition occurs not only at half fill- 
ing but also other filling cases, and in continuum model. 
This is peculiar to the attractive Fermi gases and is in 
contrast with the breakdown of the Mott insulator, which 
occurs in lattice systems only at half filling. 

This paper is organized as follows. Section II discusses 
a Hamiltonian with synthetic gauge fields and shows that 
effects of gauge fields can be incorporated into boundary 
conditions. Section IH examines Bethe ansatz equations 
and spin gap with complex spin twists. We analytically 
show that the spin gap closes at a critical value. In Sec. I 
V, we relate the real spin twist to the complex twist with 
the Landau-Dykhne method and evaluate the transition 
probability. We also comment on the ground-state de- 
cay rate. Section V summarizes the contributions and 
examines the continuum model. 



II. HAMILTONIAN 

We start by considering a Hamiltonian of two compo- 
nent atoms influenced by synthetic gauge flelds to yield 
the spin-orbit coupling. In the continuum space, the 
kinetic term becomes the covariant form: —(V + iAY , 

where A — X]a=i '•)''''' represents the synthetic 

gauge field and r° is the generator of SU(2). We then 
wish to consider the situation that the particles are 
loaded into an optical lattice. Let us assume that the 



lattice potential is so deep that the energy gap between 
first and second bands is much larger than the other ef- 
fects such as the thermal and mean-field interaction en- 
ergies per particle. If V is the lattice potential and uj is 
the Wannier function, which is localized at each lattice 
point, the hopping amplitude becomes 

= [ dxuj*{x - Ri)[-{V + iAf + V{x)]uj{x - Rj) 

dxuj*{x - R^)S-'^{x,xa)[-V^ + V{x)] 
xS{x,xq)lo{x ~ Rj) ^ S{i,j)tij, (1) 

where tij = J dxLu*{x - + V{x)]uj{x - Rj) is the 

hopping amplitude at yl = 0, 

S{i,j) = [l + iA{R,,t)-Ay\[l + iA{R, + Ay,t)-Ay\... 

N-l 

= lim T[[l + tA{R, + nAy,t)- Ay\ 



r!=0 



Pexp 



A{y, t) ■ dy 



(2) 



is the so-called Wilson line in gauge theory, and P de- 
notes the path-ordered product 22]. Here, to obtain the 
last expression of Eq. ([Ij, we used the assumption that 
the Wannier function is tightly localized at each lattice 
point. 

Although the above discussion is independent of space 
dimensions, geometry of systems, and statistics of parti- 
cles, we hereafter concentrate on a one-dimensional ring, 
and assume that the hopping is site-independent. Then, 
the Hamiltonian of two-component Fermi gases with a 
contact attractive interaction in units of the hopping 
strength is given by 



H 



-MY. 



(3) 



We note that this Hamiltonian corresponds to that an- 
alyzed in Ref. 23] except for points that the attractive 
case is concerned and there is no U(l) charge fiux in our 
treatment. In the following, for simplicity we consider the 
case that one of Aq, say, A^ has a nonzero value, which 
is comparable to the case of S0(2) spin symmetry. This 
simplified gauge configuration is also considered in the 
continuum space in Refs. |23. [25l . Then, by considering 
the following transformation |23| : 

= [S{i, I - \)S(i - 1, ^ - 2) • • • 5(2, l)],^c,^ , (4) 

Eq. ([3]) is rewritten as follows: 

7J = - ^ cfo-Ci+io- + H.c. - 4lu| Yj nitn-ii (5) 




FIG. 1: (Color online) Schematic illustration of spin currents 
induced by the time-dependent gauge flux ^{t) — Ft on the 
spin sector. Each black dot without an arrow indicates a 
particle with spin-up and counterclockwise flow or spin-down 
and clockwise flow. 



with 



CL+ia = exp[ia^L]cic 



(6) 



where L is the number of sites, ct = ±1 (or a represents 
t or 4-), and <f> is the gauge flux per site. Equation ^ 
indicates that the effect of gauge fields is incorporated 
into twisted boundary conditions. 

Let us briefly comment on a role of <&. As discussed 
in Refs. [H, Hy, [13] , in the repulsive Hubbard model, $ 
induces spin currents since $ affects the spin sector and 
there is no spin gap. In contrast, in the attractive Hub- 
bard model without spin imbalance, since there is a spin 
gap, spin currents do not flow with a static $, and there- 
fore does not play any role. In subsequent sections, 
however, we show that the time-dependent gauge flux 
<I>(t) — Ft induces spin currents through a many-body 
quantum tunneling as shown in Fig. [TJ 



III. BETHE ANSATZ EQUATIONS WITH 
COMPLEX SPIN TWISTS 

Let us consider the Hamiltonian ([5]) with complex spin 
twists, which is equivalent to the substitution ^ — i'i! in 
Eq. ^ and is related to the repulsive Hubbard ring with 
complex charge twists. As for complex charge twists, 
this model has been phenomenologically introduced in 
Ref. [3| to describe the dissipative tunneling into the 
environment and analyze the breakdown of the Mott in- 
sulator. In fact, it has been shown that due to complex 
charge twists the Mott gap closes, which evidences the 
breakdown of the Mott insulator. In this section, we dis- 
cuss the spin-depairing transition due to complex spin 
twists. The connection between the Hubbard ring with 
real twists and that with complex twists is shown in Sec. 
IV based on the analysis in Ref. (T^ . 

Now we examine solutions of the Hubbard ring with 
complex spin twists. With absence of the twists, our 
model reduces to the one-dimensional Hubbard model 
with periodic boundary conditions, which can be solved 
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exactly by the Bethe ansatz method fi\. While in the 
presence of the boundary condition is altered from 
periodic to twisted ones, we can still obtain exact solu- 
tions by the same method flE,^, and the corresponding 
Bethe ansatz equations are then for charge kj and spin 
Aq rapidities, 



N 

n 



An. — sin ki 



N 

n 

13=1 



T A^ — sin kj + i\u\ 
Xp ~ sin kj 



lU 



(7) 



Aq, — sin kj — i\u\ 



n 

13=1 



Aq — A^ + 2i\u\ 
Aq — Ah — 2i\u\ ' 



(8) 



where A'^ and are the numbers of total particles and 
particles with down spin, respectively. We note that the 
term e^*^ in Eq. ([8]) is peculiar to the case of the spin 
twist. The energy of the system is then given by 



E = — 2y~^cosfcj. 



(9) 



In this paper, we only consider the case of = SiVj^, 
which implies that there is no spin imbalance. When 
"I/ = 0, the ground state is then filled with the spin pairs, 
each of which forms the bound state of a particle with 
spin-up and a particle with spin-down to be called the 
fc — A string [l3| : 



sin/c± = AQ±^|u| + 0(e-''^), 



(10) 



where rj is assumed to be a positive number of the order 
of unity. We see that Eq. ([T0| is consistent with Eq. ([8|) . 
When ^I^ ^ 0, however, the k — X string solutions (|10p no 
longer satisfy Eq. (jSj. At the same time, if VE* is less than 
a critical value, which is determined below, it is expected 
that k — X string solutions still exist in the ground state. 
In fact, we find that the following fc — A string trivially 
satisfies Eq. ([8]): 



sin fc^ = Aq ± i\u\ + 0{e 



-('?T*)i 



(11) 



where rj ± 'i/ are assumed to be positive. Namely, the 
modification only occurs in the order of e~^^. Taking 
into account the fact that it is difficult to show the math- 
ematical existence of Eq. (flUl) . the same difficulty also 
occurs in Eq. pT|) . However, since we can check the 
consistency with the Hubbard ring with complex charge 
twists at half filling where there is the SO (4) symme- 
try we believe that Eq. ((TT]) is correct at arbitrary 
filling. On the other hand, Eq. ([7]) becomes 



i{k++k-)L _ 



n 

p=i 



Aq — A^ — 2i\u\ 
Aq - A/3 -I- 2i\u\ 



(12) 



By taking logarithm of Eq. (IT2]) , we obtain 



N, 



2LRe [sin-i(AQ + i\u\)] = 2t:J^ - ^ 

p=i 



Aq — Af 



(13) 



where 9{x) = — 2tan~"'^ (a;/|u|) and Jq is integer or half 
integer. We note that since Eq. ([T^ does not depend 
on ^P, the configuration of spin rapidities at ^ ^ is 
equal to that at = 0. Considering that the string 
solutions require the tight condition between fc^ and Aq 
and therefore Aq cannot be changed by 5*, this may be 
the natural result. 

Let us next consider the spin excitation above the 
ground state. As stated above, the spin excitation has 
a gap to be attributed to the appearance of fermionic 
superfluidity. In order to see what happens by complex 
twists, let us consider a spin-triplet excitation, which is 
accomplished by the manipulation breaking one of the 
fc — A string pairs and creating two charge rapidities fci 
and fc2 [ISl . We find that the corresponding Bethe ansatz 
equations are given by 

2LRe [sin-i(AQ -I- i\u\)\ = 2t:J^ " I] ^ ( ^ 



/9=1 
2 

■^6I(Aq - sinfcj). 



(14) 



N<-1 



kjL = 2TrIj + iL^ - ^ 6'(sinfcj - A/j), (j = 1 or 2) 

13=1 

(15) 

where Jq and Ij are integers or half integers. Compared 
with the Bethe ansatz equations in the ground state, it 
is clear that solutions of the above equations depend on 
^ explicitly. Since Aq represent the fc — A strings, it is 
expected that the shift of the distribution of Aq by ^ 
is tiny and major contributions by come from kj. In 
fact, by solving the above equations for a finite system 
numerically, we confirmed them. This is the same sit- 
uation as in the repulsive Hubbard ring with complex 
charge twists where the behaviors with respect to kj and 
Aq are interpreted as the so-called charge-spin separation 
[l6j . As for Ini(fcj), regardless of filling this takes a max- 
imum value at Re(fcj) = as can be seen from Fig. [2] 
This contrasts to the repulsive Hubbard ring with com- 
plex charge twists where Im(fcj) takes a maximum value 
at Re(fcj) — ±7r. This difference between repulsion and 
attraction is related to the observation that at least for 
^' = 0, the charge excitation of the repulsive Hubbard 
ring corresponds to the spin excitation of the attractive 
one by the substitution fc ^ tt — fc [l^. The result of Fig. 
m implies that this still holds for 7^ 0. 



Analysis in the Bulk Limit 

We hereafter consider the bulk limit, namely, N ^ 00, 
L — >■ 00, and N/L to be a constant, and introduce the 
ground-state distribution of Aq as (To(Aq) = l/i(AQ-|_i — 
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(a)L = 2iV^ = 50 



(a)|M| = 0.25 




(b)L = 50, = 15 





(b)|M| = 0.5 




1.0 and * = 0.2. By 
comparing (a) with (b), we see that for fixed ^, the maximum 
value of Im(fc-, ) at half filling is larger than that below half 
filling. 



Aq). By taking derivative of Eq. (fT3|) . we have 



{c)\u\ = 1.0 



-Re 



2\u\ 



(16) 



This equation determines the ground-state distribution. 
We note that B is related to filling n as follows: 



- = 2 / dACTo(A). 

1^ J-B 



(17) 



It is straightforward to check that when B = oo, n = 
1, that is, half filling. Meanwhile, when B < oo, we 
obtain n < 1. In addition, by performing the particle- 
hole transformation, the results of n > 1 can be obtained 
from those of n < 1 [l^. Therefore, in what follows, we 
only analyze the cases of n < 1. 

As for the spin excitation, by introducing the counting 
functions Zc{kj) = Ij/L and Zs{Xa) = Ja/L, Eqs. 




(d)|it| = 2.0 




FIG. 3: ^Pcr in several values of |u| as a function of filling n. 
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and H^TEl reduce to 

1 1 ^ 

Zg{X) = -Re [sin~^(A + i\u\)] + yT^^^-^ ^ 



(a)|M| = 0.25 



Zc{k) 







B 








dX'6 








-B 




k 








2^ 


" 2^ 


+ 2^ J 





A - A' 



a(A'), 



(18) 



dA6l(sinfc - A)(t(A), (19) 



where a{\) = ao{X) + [a'[' {X) + a^^X)] / L with the shift 
of the distribution crJ^^(A). By taking derivative of Eq. 
(fT8| . we obtain 



a{X) 



-Re 



— F 



ttL^^ |w|2 + (A-sin%)2 



1 



^ ,y 2|»|a(A0 
TTJ-B (21^.1)2 + (A- A')2- 



(20) 



We note that this equation tells us when the analytic 
properties of the system change. In fact, there is a possi- 
bility that the last term of the right-hand side has poles 
in the complex fc plane. Since the poles of the nearest 
to the real axis are ±isinh~^ the critical value 6cr 
(= Im(A:cr)) is 



sinh 



(21) 



which is equal to the condition that the Mott insulator is 
broken in the repulsive Hubbard ring ,16]. If 6 = Im(A;) is 
less than 6cr, the system is still in the spin-gapped phase 
in which the distribution is obtained from Eq. ([T6|. In 
contrast, at the critical value 5cr, the analytic properties 
of the system change. Therefore the breakdown of the 
spin-gapped phase, namely, the spin-depairing transition 
occurs. 

A relation between ^ and h is obtained from Eq. ([TS]) . 
By using the fact that the maximum of Im(A:) is located 
on Re(fc) — and zjyk = ib) = 0, and neglecting 1/L 
corrections on the distribution, we have 



dA6l(A + i sinh5)cro(A), 



(22) 



which determines the critical value of the twist, ^I^cr- We 
note that the above equation is equal to that appearing 
in the breakdown of the Mott insulator in the repulsive 
case [iB] if -B = oo, which corresponds to the case without 
imbalance. Figure [3] shows in several values of |u| as 
a function of filling and represents that VPcr increase with 
decreasing filling in any This implies that the spin- 
gapped state becomes robust as filling is decreased, and 
is consistent with Fig. [21 which indicates that b decreases 
with decreasing filling. In addition, while the increase of 




(b)|M| = 0.5 




(c)|it| = 1.0 




(d)|?i| = 2.0 




FIG. 4: Spin gap as a function of filling n at ^' = 0. 



^cr for small |u| is tiny at large filling and sharply-rising 
at small filling, that for large |m| is always smooth. This 
is because for large the main contribution of Eq. ([2^ 
comes from the first term but for small \u\ that is not so. 
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Spin Gap 



Gap Closing at the Critical Twist 



We now evaluate the spin gap based on the above ar- 
guments. Ahhough the analytical expressions can be ob- 
tained at half filling with the Fourier transformation as in 
Ref. (l6l |. we wish to consider arbitrary filling cases and 
continuum model. To this end, we adopt an approach 
with the dressed energy [l^ [2l|, where the elementary 
excitations over the ground state are obtained by solving 
the following equations: 



e(A) = -4RcVl-(A-i|it|)2 - 2^ 
f /-^ .., 2\u\e{\') 



dX' 



TT J-B 



(2|u|)2 + (A-A')2' 



(23) 



e(fc) — —2 cos k — fi 



1 



dX- 



\u\eiX) 



(sinfc- A)2 



(24) 



Here, /i is the chemical potential, — e(A) denotes the 
charge (pair) excitation, which satisfies e(±i3) = 0, and 
e{k) denotes the spin (particle) excitation. We note that 
when k — ibcr, Eq. (|24p is nonanalytic, which is consis- 
tent with the argument below Eq. (|20| . This is due to 
the fact that Eqs. (gSI) and ^ are related to Eqs. ^ 
and (Hi]) [U. At hah fihing, /i = -2\u\ and B ^ oo, and 
it is straightforward to check that in this case, — e(A) and 
e{k) correspond to the expressions of the charge and spin 
excitations at half filling, respectively. When B = 0, the 
system is in the empty band, and the chemical potential 
is then given by —2-^/1 + Therefore, possible values 
of chemical potential are 



- 2^/1 + |u|2 < ^ < ~2\u\. (25) 
The spin-triplet excitation is then given by 



AE{ki, k2) = e{ki) + e{k2) = -2^ - ^ 2cosk, 

j=i - 

\u\<X) 
_B -I- (sinfcj - A)2 



(26) 



Figure HI depicts the spin gap with 5* = in several values 
of |u| as a function of n and shows that the spin-gapped 
state is robust at small n with the fixed This is 
because the relative effect of \u\ is significant and it is 
expected that the tightly-bound spin pairs are formed at 
small n. In the numerical calculation, we first determine 
e(A) by solving Eq. P5)) with the iteration method, and 
then calculate the spin gap by assigning values to Eq. 

. By substituting kj = ib into Eq. , the spin gap 
at 5* < is given by 



Aspi„(fe) 



-2/i 
2 

TT 



4 cosh b 

B 

dX— 

B 



\u\e{X) 



\u\^ + [ismhb- Xy 



As discussed above, since the analytic properties of the 
ground state change at the critical twist, the spin gap 
should close at the value. We analytically show this as 
follows. When 6 = 6cr — 0, we have 



|up -f (isinh6 - A)^ 



2 P-"-A 



7r^(A) 



2\u\-ri + iX ' 



(28) 



where < <C 1, p. v. denotes the principal value, and 
we used l/(x -I- iij) = p.v.(l/x) — iTTS{x). By taking into 
account the property that the integral including the prin- 
cipal value vanishes, we obtain 



^spm 



(5cr - 0) = -2^ - 4v/l + (|u| - 77)2 - e(0) 



6(A) 



dX- , , 
g 2\u\ - rj + iX 



(29) 



On the other hand, Eq. indicates 

i-B 



e(0) = -2/1-4^/1 + |w|2 



6(A) 



B^'^2lu| +iX' 



(30) 



From Eqs. (gH) and we conclude that Aspi„(6bi ) = 
at arbitrary filling, which ensures that the spin-depairing 
transition occurs at the critical twist. 



IV. TRANSITION PROBABILITY 

As first discussed in Ref. the Landau-Dykhne 

method [l8l - l20l | makes it possible to relate the original 
Hermitian problems ([5]) and ^ to those of the Hubbard 
ring with complex twists analyzed in Sec. IH. To this end, 
let us consider a time-dependent gauge flux $(i) = Ft 
with a constant F in Eqs. ^ and ([5]), which is switched 
on at i = 0. Then, based on the analytic continuation 
of the energy levels as a function of complex time, the 
Landau-Dykhne method enables us to evaluate the spin- 
depairing transition probability, which is given by 



exp 



2 Im / di'Aspin($(t')) 







(31) 



where t* is the complex time in which the level crossing 
occurs. As pointed out in Ref. 
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the above integral 
path lies on the imaginary axis in the bulk limit. This 
implies the relationship $ = i^" in the same limit, which 
(27) reduces to the Hubbard ring with complex twists. There- 
fore, the transition probability becomes P = 



7 



with 



(a)|M| = 0.25 



2 




i^th = - / Aspin(6) — 

TT Jq do 


^ /'Sinh ^ u| 








71" io 


J-B 



db 



d\- 



-fi — 2 cosh 6 

TT 



B 

dX— 

-B \u\ 



2\u\ cosh baQ{X) 
|up + {isinhb- Xy 

\u\e{X) 



{isinhb — A)^ 



db, 
(32) 

where Fth is the threshold field. Figure [5] plots the be- 
havior of _Fth as a function of filling and shows that the 
system is vulnerable to the transition toward half filling. 
In addition to this, the increase of Fth for large \u\ is 
smooth, while that for small |u| has a large curvature 
around some filling. This may be because while for large 
\u\ the spin pairs are tightly bounded at arbitrary filling, 
for small \u\ the spin pairs are weakly bounded around 
half filling but tightly bounded below some filling. 

We finally comment on the ground-state decay rate F, 
which is is related to the transition probability as follows 



T/L 



aF 
2^ 



ln[l -P], 



(33) 



with an empirical factor a to be of the order of 1 and to 
describe the suppression of the tunneling. In the repul- 
sive Hubbard case at half filling, this empirical factor has 
been attributed to the pair-annihilation processes and 
determined as a function of u with the time-dependent 
density matrix renormalization group [28| . On the other 
hand, in the attractive Hubbard case at arbitrary filling, 
a is still needed to describe the particle-annihilation pro- 
cesses but has yet to be determined. The determination 
of a may be done in the same way as [28j . 



V. SUMMARY AND OUTLOOK 

We have investigated the attractive Hubbard ring with 
a time-dependent gauge flux on the spin sector combin- 
ing the methods of the Bethe ansatz with complex twists 
and Landau-Dykhne, and shown that the spin-depairing 
transition occurs. We point out that this transition may 
represent the nonequilibrium transition from fermionic 
superfluids to normal states with spin currents triggered 
by a many-body quantum tunneling. By using the for- 
mulation with the dressed energies, it has been shown 
that as the spin gap increases with decreasing flUing, the 
threshold field strength in which the spin-depairing tran- 
sition occurs becomes large as filling is decreased. Al- 
though our analysis has been performed for the specific 
time-dependent gauge flux as $(t) — Ft, taking into ac- 
count the fact that the breakdown of the Mott insulator 
occurs for another time-dependent gauge flux [ll], the 
same consideration may be obtained in the spin-depairing 
transition. 




3.0- 



2.0- 



1.5^ 
0.0 



(b)|?i| = 0.5 




{c)\u\ = 1.0 




(d)|?t| = 2.0 



0.2 



0.6 



0.8 



1.0 



FIG. 5: Value of the threshold field strength. In (a), Fth / 
in every filling while it takes a small value. 



We note that the analysis in the Hubbard ring so far 
can be directly applied to the case of a continuum model, 

H= f dx[-J2^Ux) {dec + iaAix,t)f (j),{x) 

-4|u|^;(a;)^|(x)0;(x)0^(x)j. (34) 
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This is because effects of gauge fields are again replaced 
by boundary conditions [29|. The continuum model is 
reached by passing to limits |m| — > and n <C 1 in the 
Hubbard model, which is due to the fact that mathemat- 
ically the continuum model corresponds to the contin- 
uum limit of the Hubbard model 13]. By considering the 
above limits, we have in the continuum model, sinh — > 6 
in Eq. 1^, Re^l - {X ~ i\u\)^ ^ 1 - {X^ + |wP)/2 in 
Eq. ([23)) . and cosfc — > 1 - and sinfc — > fc in Eq. 
(1211). Therefore, we do not need a special consideration 
to calculate the transition probability in the continuum 
model. 

We also point out that although we have analyzed the 
case of the spin SO (2) gauge field for experimental sim- 
plicity, the same treatment is applicable to cases of SU(2) 
gauge fields since effects of the gauge fields are again in- 
corporated into boundary conditions As far as the 
continuum models are concerned, in general, we obtain 
exact solutions by the Bethc ansatz method for attractive 
SU(A'^) Fermi gases on a ring in which the spin excita- 



tions have gaps [s^. Hence, the spin-depairing transition 
in the SU(A^) Fermi gases is also expected, which may be 
realized with, for example, ytterbium fermionic isotopes 

m. 
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